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Abstract 



Brillouin scattering studies have been carried out on high quality single crys- 
tals of Fes04 with [100] and [110] faces in the temperature range of 300 to 
30 K. The room temperature spectrum shows a surface Rayleigh wave (SRW) 
mode at 8 GHz and a longitudinal acoustic (LA) mode at 60 GHz. The SRW 
mode frequency shows a minimum at the Verwey transition temperature Ty 
of 123 K. The softening of the SRW mode frequency from about 250 K to Ty 
can be quantitatively understood as a result of a decrease in the shear elastic 
constant C44, arising from the coupling of shear strain to charge fluctuations. 
On the other hand, the LA mode frequency does not show any significant 
change around Ty, but shows a large change in its intensity. The latter shows 
a maximum at around 120 K in the cooling run and at 165 K in the heating 
run, exhibiting a large hysteresis of 45 K. This significant change in intensity 
may be related to the presence of stress-induced ordering of Fe 3+ and Fe 2+ 
at the octahedral sites, as well as to stress-induced domain wall motion. 

PACS No.: 78.35.+C, 62.65.+k, 71.30.+h 
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I. INTRODUCTION 



Magnetite has attracted renewed attention in recent years because of its interesting elec- 
tronic, magnetic, and transport properties as well as potential applications in the magnetic 
multilayer devices [1-4]. Stoichiometric Fe 3 4 shows a first order metal-insulator (MI) tran- 
sition, referred to as Verwey transition, with the reduction in the electrical conductivity of 
nearly two orders of magnitude at the transition temperature (TV). The cubic inverse spinel 
structure (F3dm) of Fe 3 4 may be written as Fe 3+ [Fe 2+ Fe 3+ ]0 4 , where the tetrahedrally 
coordinated Fe 3+ ions (A-site) is differentiated from the bracketed Fe 2+ and Fe 3+ octahedral 
sites (B-sites). The spins of the two B-site cations are antiparallel to the single A-site Fe 3+ 
cation, giving rise to ferrimagnetic order, with a transition temperature of 850 K [5] . Cubic 
magnetite at room temperature changes to a monoclinic structure below Ty, containing four 
rhombohedrally distorted cubic cells [6,7]. However, there is also an evidence of a triclinic 
structure at low temperatures [8,9]. The question of charge ordering (CO) of the Fe 2+ and 
Fe 3+ states on the B sites in the low temperature phase is a matter of continued contro- 
versy. The recent experimental data on Mossbauer [10], neutron [11,12], x-ray diffraction 
[6] and 57 Fe NMR [13] studies are not consistent with the expected CO. On the other hand, 
the high-energy transmission electron diffraction studies show charge ordering below the 
Verwey transition [7]. Also, the high resolution x-ray and neutron powder diffraction data 
has revealed charge ordering, two Fe ions with a charge of +2.4 and the other two of +2.6, 
as evidenced by the presence of two types of B0 6 octahedra and the distribution of B-B 
distances [14,15]. 

The effect of nonstoichiometry, either by oxygen deficiency or impurities has been studied 
recently [16-22]. Nonstoichiometry shifts the Verwey transition to lower temperatures as 
well as broadens the transition [16]. The nature of the transition also changes from first order 
to higher order beyond a certain doping level [2]. The influence of the linear magnetoelastic 
interaction on the ultrasonic sound velocity and attenuation in Fes0 4 was investigated by 
Moran and Liithi [23]. Of the three elastic constants, only C 44 showed a significant softening 
from room temperature to Ty, whereas Cn and C12 showed an increase as T was lowered, 
with a jump at Ty [24]. 

The role of lattice in the Verwey transition in Fea0 4 has been discussed in the literature 
[12]. Recent Raman studies have revealed a strong electron phonon coupling for the T2 9 
optical phonon [25]. The motivation of the present work is to study the behavior of acoustic 
phonons across the Verwey transition using Brillouin scattering experiments. In this paper, 
we present the temperature dependence of the surface as well as bulk acoustic phonons in 
pristine Fes0 4 from 300 to 30 K. A pronounced softening of the surface acoustic phonon 
velocity as T is lowered to Ty, and its subsequent significant hardening below the transition, 
is observed in our experiments. On the other hand, there is no appreciable change in the bulk 
acoustic phonon velocity with temperature along both [100] and [110] directions. However, 
a large change in the intensity for the bulk mode is observed near Ty, with a hysteresis of 
about 45 K. This non-monotonic behavior of the intensity is may be due to the large domain 
wall movements expected close to Ty as well as stress induced ordering of the ions. 
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II. EXPERIMENTAL DETAILS 



Single crystals of magnetite used in this study were grown by the standard skull melt- 
ing technique [26,27]. Standard four probe resistivity measurements show the MI Verwey 
transition to occur at 123 K [28]. The crystals were mechanically polished using diamond 
suspension to obtain a highly reflecting surface. After polishing, the crystals were cleaned 
in methanol, followed by acetone, using an ultrasonic bath. The measurements have been 
performed on the [100] and [110] faced crystals of the undoped Fe 3 4 . Brillouin spectra 
were recorded in the oblique-incidence as depicted [29,30] 180°-back-scattering geometry, 
with an incident angle of 9 — 45° relative to the surface normal. Ideally 9 should be 70° 
for an oblique-incidence experiment. In our experiment, the horizontal plane defines the 
scattering plane and the sample is mounted with its surface perpendicular to the scattering 
plane. The optical windows in our cryostat used for the temperature dependent studies 
are placed at 90° to each other. Thus the angle of incidence has to be 45°, in order to 
avoid hitting the reflected light inside the cryostat. The scattered light was analyzed using 
a JRS Scientific instruments (3+3)-pass tandem Fabry-Perot interferometer equipped with 
a photo avalanche diode detector. The 532 nm (A) line of Nd:YAG single mode solid-state 
diode-pump frequency doubled laser (Model DPSS 532-400, Coherent Inc. USA) was used 
as the excitation source. The temperature-dependent measurements with a temperature 
accuracy of ±1 K were carried out inside a closed cycle helium cryostat (CTI Cryogenics, 
USA). The incident laser power was kept at ~ 25 mW focused to a spot of diameter ~ 
30 fim. The typical time required per spectrum is 1.5 h. The lineshape parameters - peak 
position, full-width at half maximum (FWHM), and area were extracted by non-linear least 
square fitting of the data with Lorentzian functions, together with appropriate background 
corrections. 



III. RESULTS AND DISCUSSION 

Figure 1 shows the room temperature Brillouin spectrum of Fea04 [100], displaying two 
well resolved peaks, at 8 GHz and 60 GHz. The 8 GHz mode shows a linear dependence of 
its frequency on the wavevector component parallel to the surfaces, q\\ (= 47rsin6>/A, where 9 
is the scattering angle), expected for a surface Rayleigh wave (SRW) (as shown in inset (a) of 
Fig. 1). Further, the frequency of the 8 GHz mode shows a sinosidal dependence (see inset 
(b) of Fig. 1) on the azimuthal angle <fi (angle of rotation with respect to the scattering plane 
about the normal), confirming it to be a SRW mode. For the [100] oriented crystal, the high 
frequency mode appears in HH polarization, showing it to be a bulk longitudinal acoustic 
(LA) mode. Here H stands for horizontal polarization, which lies in the scattering plane. The 
surface Rayleigh wave velocity (obtained from the slope of the curve shown in inset (a) of Fig. 
1) is 3200 m/sec, and the bulk longitudinal acoustic velocity is 6600 m/sec, in close agreement 
with the earlier reported longitudinal sound velocity along [100] direction [23]. We point that 
a very weak mode is seen at 13 GHz in unpolarized spectrum, which can be attributed to the 
magnon. This mode could not be observed with the sample inside the cryostat. We further 
note that transverse acoustic phonons, though seen in oblique incidence Brillouin scattering 
experiments [30,31] are too week to be observed in our backscattering experiments. We 
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believe that the strong absorption at 2.32 eV excitation energy [32] broadens the Brillouin 
line making it difficult to observe. 

The temperature-dependence of the SRW frequency, v s , along [100] direction of Fe 3 4 is 
shown in Fig. 2(a) . Upon cooling below 250 K, there is a large softening in u s (~ 6 %) till 
Ty. Below Ty there is a steep increase in v s (~ 10 %). The temperature-dependence of SRW 
along [110] direction is similar to that observed in [100] direction as shown in Fig. 2(b), 
though the magnitude of softening is smaller compared to the [100] crystal. The magnitude 
of anisotropy in v s in the two directions [100] and [110] is shown as a function of temperature 
in the inset of Fig. 2(b). The full width at half maxima (FWHM) as well as intensity of 
the SRW mode for both directions [100] and [110] show changes near Ty. The FWHM of 
the surface mode is comparable to that of the instrumental broadening and hence cannot be 
qualitatively interpreted. However, we note that below Ty, the full width at half maxima 
of the SRW increases, which is still to be understood. The intensity of the SRW mode 
decreases as the temperature is lowered due to the Bose-Einstein population factor. The 
SRW mode could not be observed below 50 K in the cooling cycle, but it reappears upon 
heating the sample above 50 K. 

The bulk LA velocity does not undergo any observable change with temperature, i.e., 
there is no appreciable change in the Brillouin peak position (as shown in Fig. 3(a)) or 
FWHM. However, the peak intensity shows a strong anomaly around Ty for both the [100] as 
well as the [110] directions, as can be seen in Fig. 3(b) (for the [100] direction). In the cooling 
run, the intensity gradually increases with a decrease in temperature and shows a maximum 
at ~ 120 K, followed by a sharp drop below Ty (see filled circles in Fig. 3(b)); the mode 
cannot be observed below 70 K. In the heating run, the LA mode reappears at ~ 70 K and 
its intensity is a maximum at 165 K (see open circles in Fig. 3(b)), displaying a considerable 
hysteresis of ~ 45 K (Fig. 3(b)). The hysteresis is approximately double that which is 
measured in the magnetization measurements [33]. We now explain our observations. 

A. Surface Rayleigh wave 

For the cubic crystal, the SRW frequency depends on the elastic constants Cn, C12 and 
C44 [34,35]. Using Green's function method [36,37] we have calculated surface mode spectra 
for different temperatures in the [100] direction with the elastic constants Cn, C12 and C44 
taken from Ref. [24], and shown by dotted lines in the inset of Fig. 2(a). The values of v s 
thus calculated are plotted as dotted line in Fig. 2(a). Though the temperature dependence 
of the calculated v s is similar to our observation, the quantitative agreement is not so good. 
It is important to note that the values of the elastic constants shown in the inset of Fig. 2(a) 
were obtained in a magnetic field of 0.5 T. In order to obtain a better fit of the calculated 
v s with the data, we adjusted the value of C 44 , keeping the values of Cn and C12 as before. 
This was done because the temperature dependence of v s is similar to that of C44. The solid 
line in Fig. 2(a) shows the calculated v s , in very good agreement with the data. The value 
of C44 needed to get this good fit are shown by filled circles in the inset of Fig. 2(a). 

We will now address the softening of C 44 between 250 K and Ty. This behavior of C 44 
has been attributed to a bilinear coupling of the shear strain, e, with the order parameter, 
77, related to the charge ordering processes [23,24]. A similar behavior of the elastic constant 
C44 has been observed in other charge ordered systems Yb4As 3 [38] and a'-~NaV 2 5 [39]. 
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Taking the usual Landau free energy expansion in terms of an order parameter 77, along with 
the additional strain-order parameter coupling term Fc leads to the free energy expression 
F = F Q + arf + firf + F c + F d . Here F is the part of the free energy which remains constant 
through the phase transition, F e[ is due to the elastic energy, and F c = gerj, where g is 
the coupling constant between the strain and the charge fluctuation order parameter. This 
coupling leads to the softening of the shear elastic constant [38] 

C 44 (T)=C U (T-T )/(T-Q) (3.1) 

where 6 is the transition temperature without the strain-order parameter coupling, and 

T = 6 + 0.5#VC 4 4 (3.2) 

with the usual Landau term a = a'(T — 0). The difference of the two characteristic temper- 
atures (7q — 6) corresponds to the coupling strength of the strain and the order parameter 
[38]. The values of C 44 , obtained by the fit of the SRW frequency at various temperatures 
shown by the filled circle in the inset of Fig. 2(a) were fitted to give T = 57 K and = 
40 K, showing the coupling strength between strain and order parameter to be around 17 
K, which is somewhat higher than that estimated (~ 10 K) by Schwenk et al [24]. Since 
the contribution of C44 to v s in [100] and [110] directions are different, the anisotropy pa- 
rameter plotted in the inset in Fig. 2(b), will also reflect the temperature dependence of 
C44. The gradual decrease in intensity with temperature is expected from the usual thermal 
population factor. 

B. Bulk LA mode 

Unlike transparent materials, absorbing materials have a complex refractive index, 
n = r\ + in and the optical field falls off exponetially in the medium. This results in 
the uncertainity in the wavevector transfer q [40-42] , leading to broadening of the Brillouin 
lines, over and above the lifetime broadening. The full-width at half maximum (FWHM), 
T = AvkoK and the peak position at uj = 2r/vk gives £ = 2^; where v is the sound velocity, 
k Q is the incident wavevector [29,43]. For Fes0 4 , a has been measured only up to 0.75 eV 
and its value is ~ 3.5 x 10 3 cm" 1 [32]. Unfortunately, there is no measurement of a at 2.32 
eV. Extrapolating the results of Balberg et al [32], we find a ~ 10 5 cm -1 at 2.32 eV. This 
gives T = 14 GHz, which is close to the observed value of 16 GHz in our experiments. Hence, 
most of the width of the bulk LA mode is due to the large absorption coefficient of Fes04. 

Schwenk et al. [24] have shown that the change in Cn is ~3 % between 300 to 80 K 
(see inset in Fig. 2(a)), which implies a change of 1.5 % (~ 0.9 GHz) in LA frequency, 
since the change in density is negligible [44] . This is within the resolution limit of 1 GHz in 
our Brillouin experiments and is consistent with the temperature dependence shown in Fig. 
3(a). Now we give a plausible explaination of the non-monotonic variation of intensity of 
the LA mode with temperature (Fig. 3(b)). The intensity of the scattered light by the LA 
mode is related to the Pockels coefficients pn and p± 2 , which are related to the derivatives 
of the dielectric function with respect to the strain e. In the case of Fe 3 4 , the strains are 
not completely elastic. There are two sources of non-elastic strains; one is a stress-induced 
domain wall motion (magnetostriction) [45,46], and the other is stress-induced ordering 



5 



of Fe 2+ and Fe 3+ among the octahedral sites, e rder- Hence the total strain e is the sum 
of ^elastic, Q and border- Fine et al. [47] have found that the Young's modulus Y (Y — 
stress/e) of magnetite shows a minimum near Ty (see inset of Fig. 3(b)), which is absent in 
a magnetic field. Since a magnetic field should have no effect on elastic strains and should 
not prevent ordering, the anomaly in the Young's modulus is attributed to the contribution 
of the domain wall motion to the strain. Bickford et al. [48] suggested that if a single crystal 
of Fe3C>4 is cooled through the transition in absence of magnetic field, it is not a true single 
crystal below Ty . Furthermore, the regions may be twinned, since the a and b crystal axes 
can be interchanged. In some of the earlier experiments the domain size has been found 
to be around 0.5 up to 2 fim [49,50]. This is much smaller compared to the laser spot of 
30 /im diameter. However, the crystal reverts to the single crystalline cubic phase upon 
heating above Ty [48]. An increase in strains around Ty is expected to result in an increase 
in the value of the Pockel coefficients and hence, the intensity of the LA mode. It would 
be interesting to measure the intensity of the mode in the presence of magnetic field, which 
will make the sample single domain. 

In conclusion, Brillouin scattering studies on Fe3C>4 reveal that there is a large change in 
the surface Rayleigh wave velocity associated with the Verwey transition, which is driven by 
a bilinear coupling of shear strains with the charge ordering. The transition temperatures 
T and 9 are found to be 57 K and 40 K, showing the coupling strength to be 17 K. The 
temperature behavior of the intensities of the longitudinal acoustic mode may be understood 
in terms of non-elastic strains due to domain wall motion and charge ordering. 
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FIGURES 



FIG. 1. Room temperature Brillouin spectrum of Fes04 along [100] direction showing both the 
surface Rayleigh wave (SRW) at 8 GHz and bulk mode (B) at 55 GHz. Inset (a) and (b) shows 
the wave vector dependence and angular variation of SRW mode velocity, respectively. 

FIG. 2. Temperature dependence of the SRW mode frequency (a) along [100] direction and 
(b) along [110] direction for Fes04. The dotted line in (a) is the calculated Brillouin shifts using 
the Green's function, taking Cn, C12 and C44 from Ref. [24] shown by dotted lines in the inset. 
The solid line in (a) is calculated as a best fit to the data in Ref. [24], but varying C44 (as shown 
by the filled circle in the inset). The solid line in the inset is a fit to C44 (filled circle) using Eq. 
(1). Inset in (b) shows the anisotropy in frequency along [100] and [110] directions as a function 
of temperature. 

FIG. 3. Temperature dependence of the bulk LA mode (a) frequency and (b) intensity of Fe304 
along the [100] direction. The solid (open) circles represent cooling (heating) run data. The solid 
lines drawn in (b) are a guide to the eye. The inset in (b) shows the Young's modulus along [100], 
reproduced from Ref. [24]. 
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